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ABSTRACT 

The estimation of the reliability R = P (Y<X), when X and Y are independent and identically 
distributed random variables, has been extensively studied in the literature. Distribution functions 
such as Normal, Burr Type X, Exponential, Gamma, Weibull, Logistic, and Extreme Value have been 
considered. The expressions to estimate R have been expressed as function of positive parameters 
from the assumed distributions. No article yet has considered the expression for R as a random 
variable by itself. In this paper we will do so and consider some distributions that have been 
introduced in the literature to estimate it. The objective of this paper is to derive the probability 
density function (pdf) of R, under the assumption that the parameter of the distribution under 
investigation is itself a random variable with some specified pdf Once the pdf of R is derived, its 
moments and other properties are investigated. Several distributions including the Exponential, Chi- 
Square, Gamma, Burr Type X and Weibull are assumed for this research and are investigated. 
Furthermore, two cases are considered for each distribution. It has been found that the moment 
generating function techniques for finding the moments do not work. 



I. INTRODUCTION AND BACKGROUND 

The product reliability seems, finally, becoming the top priority for the third millennium and a technically 
sophisticated customer. Manufacturers and all other producing entities are sharpening their tools to satisfy such 
customer. Estimation of the reliability has become a concern for many quality control professionals and 
statisticians. Let Y represent the random value of a stress that a device (or a component) will be subjected to in 
service and X represent the strength that varies from product to product in the population of devices. The device 
fails at the instant that the stress applied to it exceeds the strength and functions successfully whenever Y is less 
than X. Then, the reliability (or the measure of reliability) R is defined as P{Y < X), that is, the probability that a 
randomly selected device functions successfully. Algebraic forms of R, for different distributions, have been 
studied in the literature. Among those distributions considered are the normal, exponential, gamma, Weibull, 
Pareto, and recently extreme value family by Nadarajah (2003). Haghighi and Shayib (2009) have considered 
the logistic and extreme value distributions. In all cases, X and Y are assumed to be independent random 
variables. If Hy(y) and fx{%) are the cumulative distribution function (cdfj and the probability density 
function (pdf) of Y and X, respectively, then it is well known that 



R = P(Y<X)=f H Y (z)f x (z)dz 



(1) 



The estimation of R is a very common concern in statistical literature in both distribution-free and 
parametric cases. Enis and Geisser (1971) studied Bayesian approach to estimate R. Different distributions have 
been assumed for the random variables X and Y. Downton (1973) and Church and Harris (1970) have discussed 
the estimation of R in the normal and gamma cases, respectively. Studies of stress-strength model and its 
generalizations have been gathered in Kotz et al. (2003). Gupta and Lvin (2005) studied the monotonicity of the 
failure rate and mean residual life function of the generalized log -normal distribution. 

Out of 12 different forms of cdf Burr (1942) introduced for modeling survival data, Burr Type X and 
Burr Type XII distributions have received an extensive attention. Awad and Gharraf (1986) used the Burr Type 
X model to simulate a comparison of three estimates for R, namely, the minimum variance unbiased, the 
maximum likelihood, and the Bayes estimators. They also studied the sensitivity of the Bayes estimator to the 
prior parameters. Ahmad et al. (1997) used the same assumptions except they further assumed that the scale 
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parameters are known. Then, they used maximum likelihood (MLE), Bayes and empirical techniques to deal 
with the estimation of R in such a case. They compared the three methods of estimation using the Monte-Carlo 
simulation. Additionally, they presented a comparison among the three estimators and some characterization of 
the distribution. The first characteristic they studied was based on the recurrence relationship between two 
successively conditional moments of a certain function of the random variable, whereas the second was given by 
the conditional variance of the same function. Surles and Padgett (1998) introduced the scaled Burr Type X 
distribution and name it as Scaled Burr Type X or the Generalized Rayleigh Distribution (GRD). They 
considered the inference on R when X and Y are independently distributed Burr Type X random variables, 
discussed the existing and new results on the estimation of R, using MLE, and introduced a significance test for 
R. They also presented Bayesian inference on R when the parameters are assumed to have independent gamma 
distribution. They, further, offered an algorithm for finding the highest posterior density interval for R and using 
different methods of estimation, calculated the value of R under the assumed distributions. 

Baklizi and Abu-Dayyeh (2006) have considered the problem of estimating R = (Y < X) when X and 
Y have independent exponential distributions with parameters and A, respectively. Also, Baklizi and El-Masri 
(2004) had considered the same problem as in Baklizi and Abu Dayyeh (2006), but with a common location 
parameter /(. Surles and Padgett (1998) considered the inference on R when X and Y are independent and 
identically distributed as Burr-type-X random variables with parameters and A, respectively. In addition, 
Raqab and Kundu (2009) considered the scaled Burr Type X distribution for comparing different estimators of 
R. Moreover, the expression for R came out as is specified in 

R = P(Y < X) = -±j, (2) 

where the parameter A is of the stress Y and is that of the strength X. It is to noted here that R increases as the 
ratio Xjd decreases, and, thus, the parameter A of the stress Y needs to be much smaller than of the strength X. 
Shayib (2005) had considered the effect of the sample size and the ratio between the parameters on the 
estimation of R, when X and Y shared the Burr Type X distribution. Haghighi sand Shayib (2008), considered 
the Weibull distribution when estimating R. A common application of the Weibull distribution is to model the 
lifetimes of components such as bearings, ceramics, capacitors, and dielectrics. Haghighi sand Shayib (2009) 
also considered the estimation of R in the case of the logistic distribution. 

Generally, the cdf and pdf of a two-parameter logistic distribution are given, respectively, as: 

x — a 




1 e 

F x (x) = ^r^ and f x {x) = - . r _ „ „ , -oo < x, a < oo, 8 > 0, (3) 



where a and Q are the location and scale parameters, respectively. We consider the unknown-one-parameter 
logistic distributions for both X and Y, that is, we assume that the location parameter a is 0. Hence, cdf and pdf 
are, respectively, as follow: 

_ v 

H Y (y) = y- and h Y (y) = y— — . -oc < y < oo, ft > 0, (4) 
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F x (x) = Tf- and f x {x) = x~^> ~°° < x <: °°> ft > °" ( 5 ) 

l+e ft 

Haghighi sand Shayib (2009) proved the following Theorem: 
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Theorem: When X and Y are independent random variables with logistic probability distribution and the 
location parameter a = 0, then 

,i=Q m=0 1 -| — 

n + lft 
They also considered the estimation of /? in the case of the extreme value distribution. They assumed 
that both X and Y are independent and each has an unknown-one-parameter extreme-value type 1 (or double 
exponential or Gumbel-type) distribution. Without loss of generality, it can be assumed that the location 
parameter is zero for each of X and Y. Hence, in this case, the cdf and pdf for Y and X are, respectively, as 
follows: 

_y_ _y_ — 

F Y (y.0 l )=e A. and/(j;^) = z-« & e~ e A , 00 < y < 00, d x > 0, (7) 



and 






F x (x; a ) = e" e ^ 2 and h(x; 2 ) = ^e e 2 e _e $2 , -00 < .3; < oo, <9 2 > 0- (8) 

Moreover, Haghighi sand Shayib (2009) proved the following Theorem: 
Theorem: Given (7) and (8), if 

fcrp < 1, A; = 0,1, 2,..- , 
then, 

II. R AS A RANDOM VARIABLE 

We have seen that /? has been expressed as a ratio of one positive parameter over the sum of two 
positive parameters. This ratio has appeared under the assumption of exponential, Burr Type X, gamma, and 
Wiebull distributions. In the cases of logistic and extreme-value distributions, an infinite series presentation of R 
has been derived; see Haghighi sand Shayib (2009). Thus, in the remaining of the sequel we will assume 

R= uh- m 

where U > 0, V > and that V and V are independent random variables. 

We will consider the following distribution functions: the exponential, Chi-Square, gamma, Burr Type 
X, and Weibull. Without loss of generality and for simplifying the presentation, the one parameter in the 
considered distributions will be denoted by 0. As a special case of the gamma distribution, we will consider the 
distributions that correspond to the exponential and Chi Square. This will be dealt with in section 3 for the 
exponential distribution as the underlying assumption for the variables that are involved in R. The Chi-Square 
distribution will be taken as the underlying one in section 4. In Section 5, we will investigate a special case of 
the gamma distribution. Conclusion remarks are summarized in the last section. 

III. EXPONENTIAL DISTRIBUTIONS FOR U AND V 

We will assume that V and V are independent and identically distributed random variables that share an 
exponential distribution, but with different parameters, namely, 

U - Expifix) and V ~ Exp(p 2 ), 

where fti and ih are positive. In this case the probability density functions are known as: 
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, , , J 4- e ft , « > 

\0, u<0, 

and 



(11) 



.reft, w > 

/ 2 ( V ) = <^ 2 ^ (12) 



0, u < 0, 

Since £/ and V are independent, their joint probability density function will be 

( it. v 
g(u,v)= <j-^e ^ /j| &> u,v>0 (13) 

0, otherwise. 

In this case,y(u, v) > 0, V(ti, v) such that u > and v > 0. We are interested in R as a function of [/ and V in 
the form (10). 

Let us define 

Ri — = h\(%v) and R 2 = u + v = h 2 {u. v). (14) 

U + V 
This choice of R yields an inverse ttansformation of the following form 

U = R 1 R 2 = hi 1 (R 1 ,R 2 )andV = R 2 (l - R t ) = h^{R u R 2 ). (15) 

The Jacobian of this transformation will be 

J - d *(- J * i\)-*- 

Thus, the joint probability density function of 7?i and i?2 is given by 




(17) 



Aflifla+ftgaCl-fli) -^ i2 1|J R 2 >0 )J R 2 (l- J R 1 )>0, 

otherwise. 
From (17), the marginal density functions for R^ and H-2- For instance, for i?i we have 

h{R 1 ) = - — rg, < i?x < 1, /5i,/3 2 >0. (18) 

[(A-A)^x + A] 

Note that when the two parameters are equal h(Ri) = 1, and Ri is uniformly distributed over (0, 1). This 
situation can arise when the stress and strength, that the component is exposed to, share the same density 
function that follow the exponential, Burr Type X, or the Weibull distribution. This in turn shows the result that 
was calculated by Shayib (2005) that the value of R is 1/2, and the first moment will not exceed that. Moreover, 
when the two parameters are equal, then R-i will have a gamma density with parameters a = 2 and ft, the 
common parameter of the assumed distributions. 

Now after establishing the probability density function for R, we will find its moments. Relation (18) is 
the density function of R. Thus, we have 
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E{B) = fi = / rh(r)dr 
'o 



f rh(r)dr = Pt0 2 [ 
Jo Jo 



o [(ft-ft)r + ft] 2 



dr 



Aft 



(ft -A) 2 



In [ — e ^2 



A 



— e "2 > 1. 

Pi 



From (19) we will have 



[E(R)} 2 = /i 2 = 



2 fl2 



## 



and from (18) we have 



(ft -A) 4 

E{R 2 )= i r 2 /l(r)dr = ftft / 
Jo Jo 



ft _&£& 

In 1 — - e ' 3 2 

ft 



[(ft-ft)r + ft] 2 



dr. 



(19) 



(20) 



(21) 



Thus, the variance of R can be found from (20) and (21). 



Subsequently, other moments can be found. It is to be noted here that the Moment Generating Function 
technique to find the moments is a forbidden task. 

IV. CHI-SQUARE DISTRIBUTION FOR BOTH V AND V 

In this section, we will assume that U and V are iid random variables that share the Chi square 
distribution when the parameters are different, i.e., with degrees of freedoms r/i and d-z, d\ ^ d,2. It can be 
shown that the pdf of R will have a Beta density function with parameters and d 2 /2,c?i/2, that is 



g(r) = < 



v fck ch 

(1 " r) W| 



o, 



< r < 1, 



otherwise;. 



(22) 



From (22), the mean and the variance are given by 



fi = 



«/, 



and 



respectively. 



a 2 = 



d L + d 2 

2di d 2 



(d 1 + d 2 + 2) {d 1 + d 2 ) 



2 • 



(23) 



(24) 



When the two Chi square distributions share the same degrees of freedom, then R will have a Beta 
distribution with parameters d/2 and d/2, with the following probability density function: 



g(r) = < 



_0, otherwise. 

In this case, the mean and variance of R are well known and are given by 



(25) 



1 a 2 ! 

jj, — - and o = — -, 

2 4rf + 4 



see Hogg and Tanis (1988), for instance. 
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We note that the moment generating function for the Beta distribution cannot be expressed in a close 
form. Thus, to find other moments, we will appeal to integration methods. 

V. GAMMA DISTRIBUTION FOR BOTH U AND V 

Although the gamma distribution has two identifying parameters, by fixing one at a time, two more distributions 
will emerge. Using the notation that if the random variable X has a gamma distribution with parameters a and j3, 
that is, A" ~ G(a, 0). It is well known that when a = 1. G will be the exponential distribution with parameter 

&■ 

Now, in G we assume ft — 1. Thus, if U and V, each has a gamma distribution with the same 
parameter, say a, we see that the pdf for both U and V will, respectively, be 



/l(«) - 



rfe^- 



0, 



< u < oo, 
otherwise. 



(26) 



and 



h(:<>) = 



r-k"- 1 



e ", < v < oo, 

otherwise. 



(27) 



Since U and V are independent, it is clear that their joint density function will be 

f 1 . (uv) - 1 e~ iu+v \ < a. v < oo, 
[r(«)J 
0, otherwise. 



(28) 



Thus, in this case the joint density functions of /?i and B-> is given by 



„a-l r 2a-l^_ ri) a-l 



h{n,r 2 ) = 



T(a)]' 



0. 



e"'' 2 , r 2 >0, < ri < 1, 
otherwise. 



In accordance with Theorem 1, Section 2.4 of Hogg and Craig (1995), the random variables i?i and B-> 
are independent. The Marginal pdf of i?i is given by 



rp 1 (1 - r 1 ) a ' x T{2a) 



, o < n < l, 



«/!<••,]-<; [r(a)] (29) 

0. otherwise. 

This pdf is a Beta distribution with parameters a and a. Thus, the mean and variance of R are given, 
respectively, by 



1 A 2 ! 

ii — - and <7 = -, 



(30) 



When the parameters of the assumed gamma distributions are different, that is, with a\ and r>2, the pdf 
of R will take the following form: 



9(7-1) = < 



<f) r (?) 



< n < 1, 
otherwise. 



(31) 
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That is, R will have a Beta distribution with parameters cti/2 and 02/2. Hence, the mean and the variance are 
given by 

/<■ - , and a 2 = — — — — , (32) 

see Hogg and Tanis 1988. 

Note that the moment generating function for the Beta distribution cannot be expressed in a close form for 
finding other moments. 

VI. WEIBULL DISTRIBUTION FOR V AND V 

In this section we consider the case that the distributions of U and V are independent Weibull random 
variables, sharing the same known shape parameter a, but unknown and different scale parameters, namely, 
U ~ Weibull (r.v, /?i) and V ~ Weibull(a, ft). The Weibull distribution is commonly used for a model for life 
lengths and in the study of breaking strengths of materials because of the properties of its failure rate function. 
Resistors used in the construction of an aircraft guidance system have life time lengths that follow a Weibull 
distribution with a = 2 and ft — 10 (with measurements in thousands of hours). 

In this paper we consider the Weibull distribution with two parameters, namely, a and ft, where a is 
the shape parameter while ft 1 /" is the scale parameter of the distribution. The probability density function of the 
Weibull distribution that we will use based on this setup is 




/(«)= S rt7 u e ' U '<*>A >0 ' (33) 

otherwise. 



v 



f(v) = I j- v a ~ l e & , v, a, ft 2 >0, (34) 

[0, otherwise. 

, v f o 2 A lh (uv) - 1 e -^ u+i32V \ a, ft, ft, «, v > 0, 
g(u,v) = < x ' "■ (35) 

I 0, otherwise. 

Thus, from (14) - (16) we have: 

9{ri,r 2 ) 

o? ft ft n rl a ~ l (1 - r 1 )°- 1 e - r » l>if +A<i-i>- 1 ] , «, ft, ft, r 2 > 0. < n < 1, ^ 
0. otherwise. 

When /3i = /3 2 = /3, (36) becomes 

, , J « 2 /3 2 nr 2 "- 1 (1 - nr^er^ [t+f 1 -^)"- 1 ] . a , ft r 2 > 0, < n < 1, /<W 

f/('/-i,'''2) = \ 2 ^ ' K (37) 

0, otherwise. 



Hence, 



, . . ar a ~ l (l - r) a ~ l 

hlr) — ? . < r < 1. (38) 

[r« + (l-r; '" 



Note: 
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(i) For a = 1, certainly, R is uniformly distributed over the interval (0, 1). Hence, the mean and the variance 
of R, in this case, are, respectively, given by 1/2 and 1/12. 



(ii) Whena = 2, we obtain 



2r(l-r) 



I r ,*, 0<r < 1, 

0, otherwise. 



(39) 



It can be shown that h(r) given in (39) is a pdf f or < r < 1 . Hence, the mean and the variance of R are 
given as 1/2 and (?r — 3)/4, respectively. 

(iii) In general, for all values of a, we have the following Theorem: 



Theorem: The function in (38), for anyo , is a pdf over the range of r. 

Proof: 

It is clear that the function in (38) is positive for all assumed values of r. To prove that the integral of />(*'), over 
the range of r is 1, let r = sin t, < t < tt/2. Then, dr — 2 sin t cos t dt. Hence, we have the following: 

" /2 , . . , r 1 ' 2 2a(8ini) 2(Q_1) (cost) 2 ^- 1 ' sint cost , 
h(r)dr= / — — g-i dt 

77/2 2a(smt) 2a - 1 (cost) 2 "- 1 



(sin 2a t + cos 2 " t) 2 
77/2 2a{sint) 2Q - 1 (cost) 2 "- 1 



cos 4 " t(tan 2a t + l) 2 
12 2«(sint) 2Q - 1 

o cos 2Qf+1 «(tan 2af t + l) 2 

sec 2 i dt 



dt 
dt 



dt 



1x12 2a(tant) 2Q - 1 



(tan 2Q t+ l) 2 



1 



1 + tan 2a t 



7T/2 



= 1. 



VII. BURR TYPE X DISTRIBUTION FOR BOTH U AND V 

We now assume that U and V are iid random variables that share a Burr Type X distribution with a common 
parameter fl. Thus, the density functions of U and V are, respectively, as: 



/l(«) - 



'2Pue~ u \l - er u2 ) -\ u, /3 > 0, 



0. 



otherwise, 



(40) 



and 



f 2 (v) - 



2^-" (1 - e~ v f- 1 , u, /3>0, 
0, otherwise. 



(41) 



Hence, the joint density is, 
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20,9-1 



g{u,v) 



_ J 4(3 2 uve- (u+l ^ [(1 - e~ u )(1 - e"* )] , «, u, /3 > 0, 



0, otherwise, 

Once again, note that from (14) and (15), there are many values of (u, v) that will give the same value for Hi. 



(42) 



From (16) we have the joint probability density functions of /?i and Hi as 



9(ri,r 2 ) = 



4/3 2 r 1 r|(l-r 1 )e-N'-' +r i( 1 -'' 1 ) 2 ] 



(l- e - r ? r 2)(l _ e -^(l-^) a ) 



/?-l 



0. 



< n < 1, r 2 > 0, 

otherwise. 



(43) 



For /3 — 1, (43) becomes 



9{ri,r 2 ) 



\r 



xrfCl - r 1 ) e -[ r " ri+r ^ (1 - ri)2 ] , < n < l,r a > 0, 



othc 



(44) 



rwise. 



Thus, 



2r(l-r) 
h(r) = { [1* + (1 - rf 
0. 



f , < n < l,r 2 > 0, 
otherwise. 



(45) 



It is easy to show that (45) is a pdf when < r <1. Hence, again E(R) = 1/2 and the variance of R is 
(tt-3)/4. 

When 6 — 2, the function in (44), will take the form 



g(ri,r 2 ) = < 



16r 



0. 



1 rf(l-r 1 )e- r 2[ r i +(1 - ri)2 



1 - e -''i r '2 - e -^(i-'-i)' 2 4. e - r 2 2 H+(i-n) 2 

< ri < l.r 2 > 0. 
otherwise. 



(46) 



Letting K = r 2 + (1 — r\ ) , then g(r\ .. r 2 ) can be expressed 



g(ri,r 2 ) 



'16^(1 - n)e- A> 2 - 16rir|(l - r t ) e -^ K+r ^ 
-16rir|(l - r^e-rMK+d-r!) 2 ] + i6 ri r|(l - n)e _2ifr 3, 

< n < 1, r 2 > 0, 

0, otherwise. 



(47) 



Thus, in this case, 

h(r) = 8r(l-r) 

where A'i = r 2 + (1 - r) 2 . 



^ + 



I^ + 



+ 



A' 2 (i^ + r 2 ) 2 [A 1 + (l-r) 2 ] 2 4/i 2 



< r < 1, (48) 
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Theorem: The function h(r) in (48) is a pdf over the range of R. 

Proof: 

It is clear that for < r < 1, h(r) given in (48) is positive. To show that the integral is 1 we proceed as follows. 
Let us write h(r) as the following sum: 

K r ) = #i( r ) + 92{r) + g-iir) + gt(r), 

where 



8r(l-r) 8r(l-r) 



9i(r) = 772 > 92 (r) = 



8r(l - r) 2r(l - r) 

53(r) = [/f 1 + (i- r)r and ^ ^ — K?— 

It is to be noted that g-j (r) = 4h .(?•), where /z(r) is as expressed in (45). Thus, the integral of g-[ (r) over the 
range of r has 4 as its value. Moreover, g*{r) = h(r), and, hence, the integral of ry 4 (r ) over the range of r is 1. 
Now for integrals of r/ 2 (*') and 53(f), let x = r — 1. Thus, using this change of variable and the partial fractions 
technique, we will have the integrals of g 2 (r ) and ry 3 (?■ ) equal to 2. This shows that the integral of h(f) is 1 and 
that completes the proof. 

VIII. CONCLUSION 

Different distributions were considered for the investigation R = P(Y < X), when R was taken as a random 
variable itself, and the parameters of the assumed distributions were taken as random variables themselves. The 
ratio of the parameters has a greater effect on the value of R than the sample size involved, especially when the 
parameters are equal. Even though the stress and the strength share the same distribution, the parameter of the 
stress has to be very small compared with that of the strength, and, thus, the ratio will contribute positively for 
the value of R. Two cases were investigated under the taken distributions when the two parameters are equal and 
when they are not. Regardless of the assumed distribution, the mean of R is 1/2 when the parameters are equal. 
The moment generating function technique, to find the moments did not work in the cases that were discussed in 
this paper. Additional work and investigation to find the moments of R is under way. It should be noted that 
when the parameters of U and V are equal, integrating their joint pdf is not easy; obviously, it will be more 
difficult when the parameters are unequal. 
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